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Abstract 

By the study of the dynamical processes related to entropy, this work aims to create a mathematical tool to identify 
magnetic clouds (MCs) in the interplanetary space using only interplanetary magnetic field (IMF) data. Used as basis 
for an analysis methodology, the spatio-temporal entropy (STE) measures the image (recurrence plots) "structured- 
ness" in both space and time domains. Initially we worked with the Huttunen et al. 2005's dataset and studied the 41 
MCs presenting a shock wave identified before the cloud. The STE values for each Bx, By, Bz IMF time series, with 
dimension and time delay equal to one, were respectively calculated. We found higher STE values in the sheaths and 
zero STE values in some of the three components in most of the MCs (30 among 41 events). In a physically consistent 
manner, data windows of 2500 magnetic records were selected as the calculation interval for the time series. As not 
all MCs have zero STE simultaneously, we created a standardization index (an entropy index, called as EI) to allow 
joining the result of the three components. With the use of EI three not known MCs were indeed identified and then 
the MVA method allowed calculating their boundaries. Thus the EI is proposed as an auxiliary tool to identify MC 
candidates based only on IMF analysis. In a promissor condition, this methodology implemented gives basis for an 
automatic MC identification procedure and surely useful for space weather purposes.. 

Keywords: space plasmas, space electrodynamics, magnetic cloud, spatio-temporal entropy 



1. Introduction 



A subset of Interplanetary Coronal Mass Ejections (ICMEs) has simple flux rope-like magnetic fields, which 
are characterized by enhanced magnetic fields that rotate slq x yly th r ough a large a ngle. Such events named mag- 
netic clouds or MCs dBurlaga et al.L 119811 : iKlein and Burlagal Il982t iGosling , 1 19901) have received considerable at- 
tention because the magnetic field configuration is described by the force-free model (e.g., Lepping et all 1990t : 
Osherovich and Burlagal 1 19971 and references therein). When the MC is moving faster than the surrounding solar 
wind (SW), the plasma and magnetic field typically accumulate in front of it, forming a disturbed region so-called 
sheath region. The rapid d ecrease in the total SW pressure with solar distance is the main driver of the flux-rope 
radial expansion (Demouli n and Das so. 2009). However, if a MC is moving slower than the surrounding solar wind 
( Klein and Burlaga 1 1 982.: iZhan g and Bur lagaill988l : lBurlaga , 1988 ). the sheath region is difficult to be detected and 
the magnetic cloud could go unnoticed by the spacecraft. The shock wave produced bythe ICMEs is spatially greater 
than the magnetic cloud, so a spacecraft can detect only the shock wave sometimes (fSchwennil2006l) . Near 1 AU MCs 
have enormous radial sizes (0.28 AU) with an average duration of 27 h, an average peak of the magnetic field strength 
of 18 nT and an average solar wind speed o f 420 Km/s (Klein and Burlaga, 1982; Lepping and Berdichevsky, 2000|). 
As has been noted by many authors (e.g., IZwickl et all 1 19831 : iRichardson and Canei 1 19951) . sometimes individual 
signatures may not be detected in all ICMEs, because they are not present or even there are data gaps. 

A practical criterion to detect magnetic clouds is to calculate the plasma beta, defined as the ratio of the thermal 
and magnetic pressure (j3 - 2fiop/B^), to find values significantly lower than one. For this calculation, magnetic 
field and plasma observations at spacecraft are required. Many times the temperature and density data on spacecraft 
have many gaps during periods in which the plasma instruments could be saturated as a result of intense particle 
fluxes (for example. Bastille Day in the ACE spacecraft). If this condition occurs, it makes impossible to calculate 
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the plasma beta, but it is still possib le to detect the magnetic cloud using magnetometers data (e.g. iHuttunen et al 

ieves-Chinchilla et al.L l2005h . In this paper, we propose a nonlinear time series analysis as an auxiliary tool 
to identify magnetic clouds using only the magnetic field data. The result of this calculation can be interpreted as an 
index related to some dynamical features that identifies the MC. 



MC events 
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MC events were investigated bv lHuttunen et al.l (120051) using the minimum variance analysis (MVA) dSonnerup and Cahill , 

1967h to determine if they have a flux-rope structure. They identified 73 MCs in the period (1997 - 2003) and 7 cloud 
candidate observed by the ACE and WIND spacecraft in solar cycle 23. The axis of a MC can have any orienta- 
tion with respect to the ecliptic plane (Bothmerand Schwenn, 1994, 1998). This axis orientation, identified by (pc, 
the azimuthal direction in the ecliptic, and 0c, the inclination relative to the ecliptic, is also calculated by the MVA 
jBothmer and Schwennlll998l) . 

To classify MCs, eight flux rope categories are often used: 

- Bipolar MCs (low inclination and flux rope-type: SWN, SEN, NES, NWS), 0c <45° and 

- Unipolar MCs (high inclination, and flux rope-type: WNE, ESW, ENW, WSE), 0c > 45 °, 

where the meanings are S for south, N North, W west and E east. 

Those authors have also included seven "cloud candidate" events for which either the fitting with MVA was not 
successful (e.g. the eigenvalue ratio < 2 or the directional change less than 30 °) or there was large values of beta 
throughout the event. In that study the criterion to define a MC were based on the smoothness of the rotation in 
the magnetic field direction confined to one plane. Additionally, they required that a MC must have the average 
values of plasma beta less than 0.5, the maximum value of the magnetic field at least 8 nT, and the duration at least 
6 h. With the last two criteria, they aimed to remove the ambiguity of identifying th e small and weak MC s. All 



selected events were investigated by analyzing 1 h magnetic field data with the MVA (ISonnerup and CahillL Il967h . 
where MCs are identified from the smooth rotation of the magnetic field vector in the plane of the maximum variance 
dKlein and BurlagaLll982h . 



Nonlinear analysis 

The concept of entropy is fundamental for the study of several branches of the physics like statistical mechanics 
and thermodynamics. Loosely interpreted, entropy is a thermodynamic quantity which describes or quantifies the 
amount of disorder in a physical system. Therefore we can generalize this concept to characterize or quantify the 
amount of information stored in more general probability distributions of a time series. This is in part what an 
information theory is concerne d with. The aforemen t ioned t heory was de v eloped during the 19 40.5 andl950i_with 
the main contributions given bylSh annon and Weaven (ll948h : lRenvil(ll959b : lKolmogorovl(ll958b : and lSinail(ll959l) . 

As the solar wind plasma generates complex fluctuations in a spacecraft-detecte d signals, those fluctuations can 
be investigated with techniques adopted from the nonlinear dynamics theory (Re mva and Unn ikrishnan. 2 O10l) . In 
principle, it is possible to study the MCs by analyzing the time series of interplanetary magnetic field (IMF). So, the 
main goal in this work is to find some nonlinear tool that helps someone identifying the MCs in the solar wind. 

Ojeda et a l. (2005 ) studied time series that represent physical parameters of the solar wind (interplanetary magnetic 
field and velocity components in the plasma flow) recorded by the WIND spacecraft with a time resolution of 1 min. 
In that work they studied 20 MCs, 17 plasmoids (events not identified as MC), and 20 time series of equivalent time 
duration of quiet solar wind. The IMF B, and solar wind components in a time interval of 48 hour before each 
MC were studied. Under MC conditions, a feature was identified that the IMF has the tendency to present lower 
spatio-temporal entropjQ (STE) values than the B^ in other cases, such as in plasmoids and during quiet solar wind. It 
seems to be a very contrasting aspect to be analyzed. Thus in this work a more detailed study of the STE in MCs is 
carried out. The behavior of the STE in time series of the IMF components, measured by the ACE spacecraft with a 
time resolution of 16 s, is explored. 



name given by Eu gene Kononov's Visual Recurrence Analysi.s (VRA) software, not to be confused with spatio-temporal entropy image (STEI) 
<MaandZhan5l200lh 
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66 2. Methodology and Dataset 



In this work, the purpose is to present a tool using the STE in a method that can identify simply and quickly the 
occurrence of a MC at about 1 AU. The main advantage of the approach proposed is that it takes into account only the 
interplanetary magnetic field data to identify the cloud. Nevertheless, this STE method does not solve the problem of 
identifying the boundaries of the clouds. To delimit the boundaries, if this is the case, the methodology presented in 
the work of lHuttunen et alJ (12005 ) is available to be used . In order to test and validate our identification methodology, 
MCs identified by other authors will be also used. Additionally new MCs will be identified by this tool unravelling 
some features hidden in those time series of data. Our methodology uses a strictly mathematical criterion from the 
nonlinear dynamics techniques to identify MCs in the solar wind. This methodology is indeed a novelty. 



2.1. The recurrence plots in Visual Recurrence Analysis software 

It is presented a summary of the ideas expressed in the Eugene Kononov's Visual Recurrence Analysis (VRA) 
softwarq3 about recurrence plots (RPs). In order to present the ideas, some figures are used to guide the description. 

Figure [T] shows a RP for a simple sine wave, using the data file just included in VRA software. In it organized 
patterns of color characteristics are shown for the periodical signal. In order to allow a comparative view, a RP of 
white noise is shown in Figure |2] with the data file also included. With a different result, an uniform distribution of 
color characteristics is noticed for the random signal. 
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Figure 1: Recurrence plot of simple sine wave, using a data file included with VRA softwai'e. In the RP organized patterns of color characteristics 
are shown for a periodic signal (A color version is available in the electronic version). 



The RP is a relatively recent technique for the qualitative assessment of time series (lEckmann et al.L 119871) . A 



technique that allow someone detects graphically hidden patterns and structural changes in data or see similarities in 
patterns across the time series under analysis. The fundamental assumption underlying the idea is that an observable 
time series (a sequence of observations) is the manifestation of some dynamic process. 

It has been proved mathematically that one can recreate a topologically equivalent p icture of the o riginal mul- 
tidimensional system behavior by using the time series of a single observable variable (iTakensL Il98lh . The basic 
idea is that the effect of all the other (unobserved) variables is already reflected in the series of the observed output. 
Furthermore, the rules that govern the behavior of the original system can be recovered from its output. 

In the RPs a one-dimensional time series from a data file is expanded into a higher-dimensional space, in which the 
dynamic of the underlying generator takes place. This is done by a technique called "delayed coordinate embedding", 
which recreates a phase space portrait of the dynamical system under study from a single (scalar) time series. To 
expand a one-dimensional signal into an M-dimensional phase space, one substitutes each observation in the original 



^VRA v4.7 [http://nonlinear. 1 lOmb.com/vri/] 
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Figure 2: Recurrence plot of White Noise, using data file included witli VRA software. In the RP an uniform distribution of color characteristics is 
shown for a random signal(A color version is available in the electronic version). 
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signal X(t) with vector (}'(/) - {x(i), x{i-d), x(i-2d), ■ ■ ■ , x(i-{m- l)d}), where ; is the time index, m is the embedding 
dimension, li is the time delay. As aresult, we have a series of vectors Y - 3'(l),3'(2),y(3), ■ • ■ ,3'(A^ - (m - l)d), where 
is the length of the original series. 

With such reconstruction it is possible to reproduce the original system states at each time where we have an ob- 
servation of that system output. Each unknown state Z(f) at time t is approximated by a vector of delayed coordinates 
Y{t) - x{t), x(t - d), x(t - 2d), ■ ■ ■ , x{t - (m - l)d. After the Euclidean distances between all vectors are calculated, 
they are mapped to colors from the pre-defined color map and are displayed as colored pixels in their corresponding 
places (see Figure [1] for example). The RP is a graphical representation of a correlation integral. The important 
distinction (and an advantage) is that the RP, unlike the correlation integrals, preserve the temporal dependence in the 
time series, in addition to the spatial dependence. 

In RPs, if the underlying signal is truly random and has no structure, the distribution of colors is uniform and does 
not have any identifiable patterns (see Figure |2] for example). There is some determinism in the signal generator, 
which can be detected by some distinctive color distribution. For example, hot colors (yellow, red, and orange) can 
be associated with small distances between the vectors, while others colors (blue, black) may be used to show large 
distances. In this printed work colors are noticed as a grey pattern (from white to black). Therefore one can visualize 
and study the motion of the system trajectories and infer some characteristics of the dynamical system that generated 
the time series. Also, the length of diagonal line segments of the same color on the RP brings an idea about the signal 
predictability. But, RP is mostly a qualitative tool. 

For random signals, the uniform (even) distribution of colors over the entire RP is expected. The more determin- 
istic the signal, the more structured the RP. So for the purpose of comparison in Figure[T]the RP of a strictly periodic 
signal can see and in Figure|2]the RP of the white noise time series. 



115 2.2. The entropy concepts in recurrence plot: a review 

116 The RP is a visual tool for the investigation of temporal recurrence in phase sp ace (iTakens 1981 ). In the literature 

117 there are some methods to calculate the entropy in the RP and the phase space ( Takensl 1981 ). The method (STE) 



presented by Eugene Kononov's VRA software is the only one provided satisfactory results when it is applied to time 
series of the IMF. A brief review on the entropy concepts is presented here. 

The STE entropy was used to measures the image "structuredness" in a bidimensional representation, i.e., both 
in "space" and time domains. It was implemented in VRA software to quantify the order in RPs. In physical terms, 
this quantity compares the distribution of distances between all pairs of vectors in the reconstructed state space with 
that of distances between different orbits evolving in time. The result is normalized and presented as a percentage 
of "maximum" entropy (randomness). When the entropy has a value of 100% it means the absence of any structure 
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whatsoever (uniform distribution of colors, pure randomness, seen in Figure |2|. On the other hand, 0% of entropy 
impHes "perfect" structure (distinct color patterns, perfect "struct uredness" and predic tability, seen in Figure[Tli. 

Recurrence is the most important feature of chaotic systems (lEckmann et allll987i) . The popularity of RPs lies in 
the fact that their structures are visually appe aling, and that they all ow the investigation of high dimensional dynamics 
by means of a s imple two-dimensional plo t (jpacchini et al. , 12009). For a better understanding and quantification of 
the recurrences, IWebber and Zbiluli (119941) have proposed a set of quantification measures, which are mainly based 
on the statistical distribution of the line structures in the RP. Re currence quantifica t ion analysis (RQA) is a nonlinea r 
technique used to quantify the information supplied by a RP (Webbe r and Zbilutl Il994i IZbilut and Webben. Il992h . 
Recurrence variables are calculated from the upper triangular area of the recurrence plot, excluding the central diag- 
onal, because the plot is symmetrical about the main diagonal. The RQA can be used as a tool for the exploration of 
bifurcation phenomena and dynamics changes also in nonstationary and short time series. The entropy (ENT) is one 
of the recurrence variables of the RQA method. It is the Shannon information entropy for the distribution probability 
of the diagonal lines. That is: 



ENT = - 



k=L, 



£ pik)l082(p(k)), 



where L„„„ is the minimum length of diagonal Unes in RP and 



Pik)^ 



number of diagonal lines of length k in RP 
number of diagonal lines in RP 



(1) 



(2) 



The ENT can be calculated using the VRA software; but it should not be confused with the STE. 



Little et al. (2007) developed a Recurrence Period Density Entropy (RPDE) method, first it requires the embedding 
of a time series in phase space, which, according to Taken's embedding theorems, can be carried out by forming 
time-delayed vectors for each value x„ in the time series. Then, around each point in the embedded phase space, a 
recurrence neighbourhood of radius e is created. All recurrences into this neighbourhood are tracked, and the time 
interval T between recurrences is recorded in a histogram. This histogram is normalized to create an est imate of the 
recur rence period density function p(T). The normalized entropy of this density is the RIDE value H„orm dLittle et al 
2007h . 



Hnorm — —ilnTmnx) ' p(t)lnp(t). 



(3) 
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The RPDE value is a scalar in the range zero to one. For purely periodic signals, //„orm = (STE=0%) whereas 
for purely uniform white noise, H„orm = 1 (STE=100%). However, estimates obtained with this technique (RPDE) 
are diff'erent from those obtained with the STE. 

IPasan et al.l (120021) report an analysis, using the tools of nonlinear dynamics and chaos theory, of the fluctuations 
in the stress determined from simulations of shear flow of Stokesian suspensions. They also computed the STE using 
VRA for the stress. The calculated values of the STE for the shear and normal stresses were nearly zero, show- 
ing perfect s t ructur e in the data. They obser ved defin i te stru c ture in the phase-space pl o t of the s t ress c omponents 
dDasan et al.L l2002h . They cited the works of lPeacockl (ll983l) : ICarr and Schwart3(ll998l) . IPeacockl (Il983 ') presented 
a two-dimensiona l analogue of Kolmogoroy -Smirnov test, useful for analysing the distribution of data in two dimen- 
sions, as is the RP. ICarr and Schwartz! ( 119981) investigated the fluctuation phenomena in plasmas that often necessitates 
the analysis of spatio-temporal signals. It was shown how such signals can be analyzed using the biorthogonal de- 
composition, which splits them into orthogonal spatial and temporal modes. Several parameters allow one to quantify 
the weight distribution in the biorthogonal decomposition. The total energy of spatio-temporal signal u(t, j) is found 
to be equal to the sum of the eigenvalues, 



(4) 



They can define the relative energy of the m structure as 



EiuY 



(5) 



162 and the entropy of the spatio-temporal signal u(t, j) is defined as 



1 



N 



H(u) = -T-^ ^ E,„(u)logE„(u). 

" m=l 



(6) 



It describes how the energy is distributed across the A^, significant structures. Signal whose energy is concentrated in 
a single structure such that A^j = 1 will have very low entropy H{u) = 0, our H(u) = 1 if the energy is distributed 
equally among the Ns significant structures. 

The results will be presented in this paper shows the usefulness of STE implemented by Eugene Kononov's 
software. Other ways to calculate the entropy are not useful as an auxiliary tool to identify the occurrence of MCs. 
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2.3. Test of stationary time series using STE 

From an intuitive point of view, a time series is said to be stationary if there is no systematic change in mean (n o 
trend) and no systematic change in variance and if strictly periodic variations ha ve been removed (IChartfieldl l2003h . 

Trend estimation is a statistical technique to aid in the interpretation of data ( ChartfieldL l2003h . When a time series 
related to measurements of a process are treated, trend estimation can be used to make and justify statements about 
tendencies in the data. Given a set of data and the desire to produce some kind of "model" of those data (model, in 
this case, meaning a function fitted through the data), there are a variety of functions that can be chosen for the fit. 
However, if there is no prior understanding of the data, then the simplest function to fit is a straight line and thus this 
is the "default". Once it has been decided to fit a straight line, there are various ways to do so, but the most usual 
choice is a least-squares fit, equivalent to minimization of the L2 norm. If there is no global trend in time series the 
angle ("trend angle") between the straight line and the positive x axis must be zero, else the angle is not zero. The 
aim of this subsection is to calculate the STE versus the trend angle. We deal with time series with data file included 
in VRA to calculate the STE of each one versus the trend angle. 

Figure [3] shows a time series plot of Lorenz data file included in VRA. We gave trends to the series through 
rotations about the origin. After that, we calculated the STE of each temporal series. The results were included in 
row 2 presented in Table [T] in which three time series with data file included in VRA version 4.7 are shown. The 
effects of the trends in every series related with STE values are quantified. 



Lorenz System 
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Figure 3: Time series plot of Lorenz data file included in VRA (case with d = 0). Series rotated about origin {6 = -0.01 rad , 6 = 0.01 rad , 
6 = 0.0175 rad ) and the three resulting series also plotted. After that, we calculate the STE of each series. 

184 

185 We follow the same idea to cause a trend in time series for another cases. Sine and White Noise data file also 

186 included in VRA. The results were included in rows 3 and 4 of Table [1] For periodic time series (sine data file) the 

187 value of STE is always zero independently from the increasing trend. If the trend in the time series increases, the STE 

188 value decreases (see row 2 and 4 in Table[T]|. 
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Table 1: STE values related to trends for three time series with data file included in VRA 4.7. 



Series 


= rad 


e = -0.01 rad 


= 0.01 rad 


e = 0.0175 rad 


STE(Lorenz) 


73% 


30% 


29% 


0% 


STE(Sine) 


0% 


0% 


0% 


0% 


STE(White Noise) 


80% 


34% 


34% 


3% 



189 Also, the first difference applied in a time series is often enough to convert series with a trend into a stationary time 

190 series. The first-order differences of time series values x\,X2,xt,, ■ ■ ■ ,xn are given by a new series y\,y2, ■ • • ,yN-\, 

191 where yff^i - x^ - Xf^^x. The operation y, - x, - - Vx, is called the first difference and V is the difference 

192 operator Sometimes the first differencing is not enough to remove the trend in mean. We then try further differencing. 

193 For example, the second-order difference is given by: 

V^y, - Vjc, - Vx,_i = X,- 2x,_i -I- x,_2 (7) 

194 Our interest is to study variations in the STE values when a first-order differences are applied on stationary time 

195 series. In the previously studied time series (Lorenz, Sine and White Noise data file included in VRA), we calculate 

196 the first-order differences. After that, we calculated the STE of each time series and compared the results with the STE 

197 of the original series, which are shown in Table |2] The STE values are similar in both of them, i.e., for transformed 

198 (first-order differences) and untransformed time series. Thus, when the time series has no trend, the non-linear delicate 
structures ai^e not destroyed. 



Table 2: STE values related to the first-order differences in time series. 



Series 


Untransformed 


first-order differences 


STE(Lorenz) 


73% 


75% 


STE(Sine) 


0% 


0% 


STE(White Noise) 


80% 


82% 



199 

200 The entropy value is small and may tend to zero in any time series with trend. The STE value for nonstationary 

201 signals is not meaningful. If there is a trend in the time series, you might want to consider removing it by differencing 

202 the original time series before calculating the STE. Keep in mind, however, that taking the first differences may destroy 

203 the delicate nonlinear structure in the time series (if there is any). Following this idea, sometimes we differentiate the 

204 series of the interplanetary magnetic field before calculating the entropy. This is done on a trial basis after calculating 

205 the entropy of the original series. The reason is to remove the trend and calculate the true value of entropy of the 

206 series, the results are shown in Figure |4] (c) and will be discussed later The relationship between trend and STE 

207 values in magnetic clouds is discussed in this work. 

208 2.4. ACE spacecraft data 

209 The plasma particle detected by ACE arrived to the Earth after approximately 30 min. The Magnetic Field Ex- 

210 periment (MAG) on board ACE consists of twin vector fluxgate magnetometers to measure interplanetary magnetic 

211 field. The datqj contains time averages of the magnetic field over time periods Is, \6s, A-min, hourly, daily and 27 

212 days (1 Bartels rotation). Magnetic field vectors are given in the RTN, GSE, and GSM coordinate systems. The Solar 

213 Wind Electron, Proton, and Alpha Monitor (SWER4M) measures the solar wind plasma electron and ion fluxes (rates 

214 of particle flow) as functions of direction and energy. The dat£0 contains time averages of solar wind parameters over 

215 time periods 64s (ion data only), 128i (electron data only), hourly (all data), daily (all data) and 27 days (1 Bartels 

216 rotation) (all data). 

217 In this paper we use data from the magnetic field components with the time resolution of I6s in GSM and GSE 

218 coordinate systems, because we are interested in studying time series of magnetic clouds from 1000 to 5500 data 



3http://www.sri.caltech.edu/ACE/ASC/level2/lvl2DATA_MAG.html 
'*http://www.sri.caltech.edu/ACE/ASC/level2/lvl2DATA_SWEPAM.html 
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records. We think that some local informations are lost in the short time series analysis. We cannot calculate the STE 
in the VRA software for time series with more than 5500 data points (explained later in Subsection l3.1b . If we need 
to detect the boundari es of a magnetic cloud with MVA method then we use data with l/i time resolution (similar to 



Huttunen et alJ (|2005|)) for MAG and SWEPAM experiments respectively. 



223 3. Results and Discussion 

224 In the next subsection, we compare STE values between time series correspondi ng to MCs and sheath s. This 



comparison is performed only for IMF data, which were data of the MCs detected by ^ Huttunen et al. I (l2005h . From 



226 the 73 MCs identified by those researchers we have dealt with just 41 MCs: just for clouds preceded by shock waves. 

227 In the other subsections, we have studied time intervals of the solar wind, designated here as solar wind windows, 

228 with the purpose of validating a method for detecting MCs. Over all the windows the STE variations were studied. 

229 We named the previous methodology of "windowed spatio-temporal entropy". The criteria to select a specific solar 
wind window are presented in the begin of each subsection. 



231 3.1. STE calculated for 41 magnetic clouds between 1998 — 2003 

232 We processed ACE data, from March 1998 to December 2003, of the three magnetic field components {B^, By, B^) 

233 in GSM, with a time resolution of 16.?. We obtained the recurrence plots for 246 time series (3*(41 MCs + 4-1 S hts) = 

234 246). We have calculated the STE value for each temporal series with embedding dimension and time delay equal 

235 to one respectively. It may be noticed that STE value changes for diff'erent embedding parameters. For example, for 

236 Lorenz attractor (Lorenz data file included in VRA), STE is near its minimum when the correct embedding is used 

237 (dimension - 3, time delay = 16 for that particular data file). This and other results suggest that STE can be used 

238 to determine the optimal embedding parameters. We selected the same embedding and time delay, equal to one, to 

239 maintain equivalence in the calculation of entropy among all series. Table[3]presents the events in chronological order 

240 The STE values for the 41 MC events, presented in Table[3] are shown in Figure|4](a). At the top, the STE of the 

241 three IMF components B^, By and B^, plotted respectively as the " o ", " + " and " x ", for the MCs. At the bottom, the 

242 same for the sheath regions. In Figure|4](b), the histogram of STE obtained for the B. component taking into account 

243 the MCs and the sheath regions are shown. 

244 In Figure |4] (a), the STE values of the 246 time series (each separated by cloud and sheath respectively) were 

245 plotted in chronological order (as appeared in column 1 in Table O earlier for MCs and bellow for sheath regions. 

246 The STE values were calculated with embedding dimension and time delay equal to one respectively. Some MCs do 

247 not have values close to zero entropy in the three components simultaneously. Then, it is possible to find components 

248 with perfect structuredness (low STE) and absence of structure (high STE) in the same magnetic cloud. 

249 If we compare the STE values between the same components for the sheath and the cloud regions (for example, 

250 B,-sheath (STE = 56%) with STE B,-cloud (STE = 0%) in the event 1), we detect that for 5/41 (3/41) cases in 

251 the Z? f (By, B~) component(s) the STE in clouds is bigger than the one in sheath respectively. We found a decrease in 

252 entropy in the clouds present in all IMF components. Someon e can notice that a clear tendency of the MCs events 



253 to present STE with lower values, close to zer o, as observed inlOieda et al.l ( 120051) and extended in this work, which 

254 is a new feature added to the usual features ( Burlaga et al. , 198 1*) established to the MCs. The STE was zero in 



255 20/41, 21/41, 25/41 MCs to B^, By, B- respectively. The three components had 5 = 0% value in 17/41 MCs and 

256 1/41 sheath at the same time. The sheath region with STE = 0% corresponds to event number 06 in Table [3] The 

257 presence of another cloud in the sheath region will be examined later (subsection l3.2.2b . 

258 Figure |4](b) shows a histogram of STE derived from (a) for the B~ component of clouds (in grey) and sheaths. 

259 We have 37/41 of MCs with STE less than 40% and 37/41 of sheaths with STE greater than 40% respectively. We 

260 did some tests with the time series to explain the results above. First, if the Gaussian noise is removed from the 

261 signal and the STE calculated, the STE value tend to decrease in less than 5% from its initial value. Second, when 

262 a trend is removed of the series through a rotation (with the angle of slope line of best fit) the STE varies, but still 

263 the three components had S TE = 0% value in 17/41 MCs and 1 /41 sheath at the same time in time series with more 

264 than 5500 data points (VRA software limitation). These results are shown in Figure |5] and will be discussed later 

265 Third, by removing the trend through the first order difference in time series (see Figure|4](c)) we still have MCs with 

266 STE^O%. 
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Table 3: Solar Wind data studied (fromlHuttunen et al.l t2003)). 



No. 


Year 


Shock 


UT 


MC, start 


UT 


MC, stop 


UT 


01 


1998 


06 January 


13: 


:19 


07 January 


03: 


:00 


08 January 


09 


:00 


02 




03 February 


13: 


:09 


04 February 


05: 


:00 


05 February 


14 


:00 


03 




04 March 


11: 


:03 


04 March 


15: 


:00 


05 March 


21 


:00 


04 




01 May 


21: 


:11 


02 May 


12: 


:00 


03 May 


17 


:00 


05 




1 3 June 


18: 


:25 


14 June 


02: 


:00 


14 June 


24 


:00 


06 




19 August 


05: 


:30 


20 August 


08: 


:00 


21 August 


18 


:00 


07 




24 September 


23: 


:15 


25 September 


08: 


:00 


26 September 


12 


:00 


08 




18 October 


19: 


:00 


19 October 


04: 


:00 


20 October 


06 


:00 


09 




08 November 


04: 


:20 


08 November 


23: 


:00 


10 November 


01 


:00 


10 




13 November 


00: 


:53 


13 November 


04: 


:00 


14 November 


06 


:00 


11 


1999 


1 8 February 


02: 


:08 


1 8 February 


14: 


:00 


19 February 


11 


:00 


12 




16 April 


10: 


:47 


16 April 


20: 


:00 


17 April 


18 


:00 


13 




08 August 


17: 


:45 


09 August 


10: 


:00 


10 August 


14 


:00 


14 


2000 


1 1 February 


23: 


:23 


12 February 


12: 


:00 


12 February 


24 


:00 


15 




20 February 


20: 


:57 


21 February 


14: 


:00 


22 February 


12 


:00 


16 




1 1 July 


11: 


:22 


1 1 July 


23: 


:00 


13 July 


02 


:00 


17 




13 July 


09: 


:11 


13 July 


15: 


:00 


13 July 


24 


:00 


18 




15 July 


14: 


:18 


15 July 


19: 


:00 


16 July 


12 


:00 


19 




28 July 


05: 


:53 


28 July 


18: 


:00 


29 July 


10 


:00 


20 




10 August 


04: 


:07 


10 August 


20: 


:00 


1 1 August 


08 


:00 


21 




1 1 August 


18: 


:19 


12 August 


05: 


:00 


13 August 


02 


:00 


22 




17 September 


17: 


:00 


17 September 


23: 


:00 


18 September 


14: 


:00 


23 




02 October 


23: 


:58 


03 October 


15: 


:00 


04 October 


14 


:00 


24 




12 October 


21: 


:36 


13 October 


17: 


:00 


14 October 


13 


:00 


25 




28 October 


09: 


:01 


28 October 


24: 


:00 


29 October 


23 


:00 


26 




06 November 


09: 


:08 


06 November 


22: 


:00 


07 November 


15 


:00 


27 


2001 


19 March 


10: 


:12 


19 March 


22: 


:00 


21 March 


23 


:00 


28 




27 March 


17: 


:02 


27 March 


22: 


:00 


28 March 


05 


:00 


29 




1 1 April 


15: 


:18 


12 April 


10: 


:00 


13 April 


06 


:00 


30 




21 April 


15: 


:06 


21 April 


23: 


:00 


22 April 


24: 


:00 


31 




28 April 


04: 


:31 


28 April 


24: 


:00 


29 April 


13 


:00 


32 




27 May 


14: 


:17 


28 May 


11: 


:00 


29 May 


06 


:00 


33 




31 October 


12: 


:53 


31 October 


22: 


:00 


02 November 


04 


:00 


34 


2002 


23 March 


10: 


:53 


24 March 


10: 


:00 


25 March 


12 


:00 


35 




17 April 


10: 


:20 


17 April 


24: 


:00 


19 April 


01 


:00 


36 




18 May 


19: 


:44 


19 May 


04: 


:00 


19 May 


22 


:00 


37 




01 August 


23: 


:10 


02 August 


06: 


:00 


02 August 


22 


:00 


38 




30 September 


07: 


:55 


30 September 


23: 


:00 


01 October 


15 


:00 


39 


2003 


20 March 


04: 


:20 


20 March 


13: 


:00 


20 March 


22 


:00 


40 




17 August 


13: 


:41 


18 August 


06: 


:00 


19 August 


11 


:00 


41 




20 November 


07: 


:27 


20 November 


11: 


:00 


21 November 


01 


:00 
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Figure 4: (a) The STE values for the 41 MC events between 1998-2003 presented in Table[3] At the top, the STE for B_„ By and B-, the three IMF 
components, plotted respectively as the " o ", " + " and " X ". At bottom, the same as to above but for the sheath regions, (b) The histogram of STE 
obtained from (a) for the B- component taking into account the cloud and the sheath regions, (c) It is similar to (a), but with the trend removed 
through the first-order dilferences in time series, (d) It is similar to (b). The histogram of STE obtained from (c) for the B- component taking into 
account the cloud and the sheath regions. 



10 



267 Figure|4](c) is similar to panel (a), but we have eliminated the trend through the first order difference in time series. 

268 In this case most of the STE values in the sheath increased to ~ 90% in the three components. Figure |4](d) shows a 

269 histogram of STE derived from (c) for the B, component of clouds and sheaths. If we eliminate the trend in the series 

270 then the entropy increases. Then, there is not a nonlinear structure in the time series, the signal is nondeterministic. 

271 Also the STE increases in the clouds, but there are still MCs with zero STE. We find a limitation to interpret this case, 

272 the length of the series. The calculation of the STE with the VRA software, version 4.7, can not be made for time 

273 series with a size bigger than 5500 points because the STE results always in zero. It seems to be a limitation of the 

274 software by some reason not explained in its tutorial. 
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Figure 5: At middle, STE values versus time series length for B- component for the 41 MCs, where the " o " and " + " symbols coiresponds to 
original and transformed (remove the Gaussian noise and trend through a rotation) time series respectively. At bottom STE values versus time 
series length for the 41 sheaths. At top the number of events versus intervals of the time series length. The histogram helps to identify overlapping 
points of MC in the middle panel. The STE is zero to the right of the vertical line due to a software limitation. 



In Figure |5] at bottom, we have plotted the STE versus the length of the time series for the B, components for a 
sheath region. The " o " and " + " symbols correspond respectively to the original and transformed (by removing the 
Gaussian noise and "trend" through a rotation) time series. Panel at middle is similar to one at bottom but for the MC 
region. The numbers to the right of the vertical line show the total of events (17) that each point represent. At top, the 
histogram helps to identify overlapping points of MC in the middle panel. The STE is zero to the right of the vertical 
line due to a software limitation. At middle panel, we have MCs (" o " symbols or original) with zero STE values to 
the left of the vertical line. When we remove the Gaussian noise and trend through a rotation the STE values (" + " 
symbols or transformed) increase and are different from zero. 

In Figure|5] related to the top and middle plots, it is shown S TE = 0% in 17/41 (events 1 - 4, 6 - 10, 13, 16, 27, 

30, 33 - 35, 40 in TableO in the MC with more than 5500 points, in the other 9/41 MCs (events 5, 12, 21, 23 - 25, 

31, 32, 37 in Table O the STE value leaves the zero value when the time series were transformed by removing the 
Gaussian noise and "trend" through a rotation. So far, we have been seeking the causes of the large amount of MCs 
with STE = 0%. During a MC, the magnetic field strength is higher than the average, the magnetic field direction 
rotates smoothly through a large angle, then the periods with MCs present more trend in the magnetic behavior than 
the periods of sheaths or "quiet" solar wind. The trend is the princ ipal cause of th e lower STE values in MC. Also, the 



290 MCs are more structured than sheath and "quiet" solar wind COieda et al.L 120051) . Then it establishes a new feature to 

291 be considered in analyses. 

292 3.2. The Windowed Spatio-Tempoml Entropy 



From the 41 MCs studied the STE decreases inside the MCs in relation to the sheath. In lOieda et al.l (I2005h we 
found MCs with lower STE values than the "quiet" solar wind and other events in the interplanetary medium not 
identified as MCs. Then we have decided to analyze solar wind intervals of 10 days where we already knew about the 
existence of MCs, as the events reported in Table [3] As a criterion we have selected data windows of 2500 records 
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297 moving forward as 200 records at head, until the end of the time series was reached. For every data segment selected 

298 the STE was calculated, allowing to analyze the evolution of the STE at every 200 * I65 in time evolution for every 

299 2500 records in the time series. 

300 We select 2500 points because this interval represents an interval of 1 1 . 11 h, and the MCs have a smooth rotation 

301 of the magnetic field vector in the order of 1 day, where the field reaches a peak and decreases. With a temporal 

302 window size of 11.11 h and a resolution of 16 s it is possible to cover the entire range of trend in the most of 

303 MCs with dimensions larger than 24 h. The STE values are calculated every 0.89 h (STE time resolution adopted) and 

304 representing ~ 8% of the size of each temporal window. So the variation of the values of STE between two continuous 

305 windows (1 1.1 1 h) must also be in the order of ~ +8%, a sufficient condition for our initial purpose. We will be able 

306 to study in detail the calculation of the STE by temporal windows for six solar wind intervals following in the studies 

307 (4 MCs represented in Tableland other 2 cases). 

308 3.2.1. 07-January-1998 MC event 

309 We decided to study in more details the STE variations taking into account some days before and after the event 

310 number 1 in the Table[3] The time series of the B B,, and IMF components of MC has zero STE value (see, Fig- 

311 ure|4](a), the first event). We find in Figure|5]that the time series are composed of 6753 data points and, consequently, 

312 we have zero entropy when the STE was calculated using the VRA software. 

313 Based on the features analyzed in the beginning of this section, our hypothesis is to find zero STE in the period of 

314 the magnetic cloud to identify it. Then, we have selected a time window for 03 - 12 January, 1998 and take IMF three 

315 components in the Geocentric Solar Magnetospheric (GSM) coordinate system recorded by the ACE spacecraft with 

316 time resolution of \6s. We selected data windows of 2500 records moving forward as 200 records to head, until the 

317 end of the time series. The numbers of windows result in 258 and we calculate the STE of each one. 

318 In Figure|6l we show the values of STE versus time for the time series of IMF B 5,, and components in the 

319 solar wind. Additionally, we investigated if ACE detected some other event during these ten days. 




Figure 6: Values of spatio-temporal entropy on January 03 - 12, 1998, as a function of the time for times series of IMF B, , By and B- components 
in the solar wind. The thick curve represents the Entropy Index (EI) calculated over the analyzed period. The shock, the start and end of the MC 
are represented by three vertical dotted lines. 
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Cane and Richardson (l2003l) summarized the occurrence of interplanetary coronal mass ejections (ICMEs) in the 
near-Earth solar wind during 1996 - 2002, corresponding to the increasing and maximum phases of solar cycle 23. 
In particular, they give a detailed list of such events based on in situ observations. They reported two ICME in this 
interval, 07 January 01 : 00 to 08 January 22 : 00 and 09 January 07 : 00 to 10 January 08 : 00. Then we have one 
MC, one interplanetary disturbance that was not classified as MC, and a "quiet" solar wind period. Some features 
related to the ICME classified as MC in January 07 were reported by ICane and RichardsonI ( l2003b : for example 
ViCME = 400km/ s is the mean solar wind speed in the ICME; V„,ax = 410^m/s is the maximum solar wind speed in 
the post disturbance region ; B - l6nT is the mean field strength; Vt = 4^0km/s is the transit speed. The MC was 
also reported bv Huttunen et al. I (l2005h but with di f ferent size (see Table. O. 

We used the dates reported bv lHuttunen et al.l (l2005b to represent the shock, the start and end of the MC. These 
dates were represented with three vertical dotted lines in Figure |6] This unipolar MC has a flux-rope type ENW 
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and the observed angular variation of magnetic field is left-handed dHuttunen et al.l l2005h . The MVA method gives 
the eigenvalue ratio /I2//I3 - 48, the angle between the first and the last magnetic field vectors x - 160°, and the 
orientation of the axis {(pc, Oc) = (21 °, 52 °). The component was northward almost during the whole passage of 
the MC, then a magnetic storm ( Dstmrn = -llnT) is cau sed by the sheath of heated and compressed solar wind plasma 



piled up in front of CME ejecta dTsurutani et all 1 19881) 



Figure |6] shows the behavior of STE in ten days of solar wind. The By, B^ components have a zero STE value 
only during the passage of the MC, approximately in the first half of it. The second mi nimum value of S TE - 10% 
corresponds to Z? component in January, 09. At this date (09 Jan 07 : 00-10Jan08 : 00) Cane and RichardsonI ( 2003 ) 
detected one ICME with Vicme = 4-5QKm/s, V„,ax - 500Km/s, B - 6nT. This result seems to be very interesting 
becaus e in ten days of solar wind data two magnetic components had zero value only within a MC. Huttunen et alJ 
( l2005h first performed a visual inspection of the data to find the candidate MCs. This is always the first step in any 
work aiming at studying magnetic clouds. Then, the calculation of STE could be used as a mathematical tool to help 
finding the candidate MCs and being a new feature of many MCs. To generalize this idea it is necessary to study 
other MCs and explain physically two questions: (1) Why STE is zero only somewhere inside the MC? (2) Why don't 
all the magnetic components have zero STE inside the MC? These questions will be answered later. As not all the 
magnetic components have zero entropy at the same time, we created a standardization index (called by us Entropy 
Index (EI)) to allow joining the results of the three variable (B^, By, B,) based on the physical process related to the 
MC geometry. The index is the result of multiplying the values of entropy of the three variables at the same time t 
normalized by 10"*: 



(S) 



The thick curve in Figure |6] is the representation of the EI calculated over the analyzed period. We could think 
of a similarity with the Dst index, which is a geomagnetic index which monitors the world wide magnetic storm 
level. Negative Dst values indicate that a magnetic storm is in progress, the more negative Dst, the more intense the 
magnetic storm. But the entropy index value can only decrease to zero somewhere inside the MC. The MCs have 
simple flux rope-like magnetic fields, characterized by enhanced magnetic fields that rotate slowly through a large 
angle. Then, the time series of IMF components have a trend and a mor e ordered dynamic behavior and a higher 
degree of correlation between its temporary neighbors (lOieda et al.l l2005b . The above mentioned behavior is found 
only in the magnetic structures of the MCs, a necessary condition for the zero entropy index value. Conclusion, using 
the entropy index we find the MC occurrences. But not the frontiers. 
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3.2.2. 20-August- 1 998 MC event 

We take this temporal windows of ten days because we already knew that there was an MC near the center, as 
reported in Table [3]number 6 event. Also, this is the only event that has zero STE in the sheath. 

Figure |7]is similar to Figure|6 l but in the period 15-25 August 1998. At this date (19 Aug 18 : 47 - 21 Aug 20 : 
00, 1998). 1Cane and Richardson! (l2003h detected one ICME classified as MC with Vjcme = 300Km/s, V,nax = 
340Km/s, B - lAnT . The three numbere d (1,2,3) vertic al dotted lines in Figure |7] represent the shock (1), the 
start (2) and end (3) of the MC2 reported bv lHuttunen et al.l ( i2005) (see event number 6 Table[3]l. The MVA method 
toMC2 gives A2IAJ. =30,;^-= 177 ° and (<;*c, 6*0) = (1 13 °, -16 °). This bipolar MC has a flux-rope type SWN and 



the observed angular variation of the magnetic field is right-handed (^ Huttunen et al. , 2005). The observed SN-type 
MC had < for about eleven hours, with values of less than -lOnT . In the sheath it was observed B, < for 
more than seven hours and values close to B, - - WnT . The Dst index shows two geomagnetic storms, the first with 
Dst„,j„ - -42nT caused by the sheath, the second with Dst„,i„ - -61nT caused by the MC. 

Figure Q shows that the B^, By components have a zero STE value during the sheath of the MC, approximately 
at the end half of it. After that, approximately in the center of the MC, the three components (Bi, By, Br) have zero 
STE values. The EI is zero at the end of the sheath and within the MC. We think that maybe there are two magnetic 
clouds, one after the other. It is very difficult to detect the two MCs when a visual in spection of the data is performed 
and then apply the method of minimum variation (see Bothmer and SchwennL 1998 ). Someone could make a mistake 
in the interpretation of the cloud considering it as part of the sheath of the other MC. To solve the above hypothesis, 
we used the same method (MVA) explained in .Huttunen et alj (l2005h . with success. 
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Figure 7: Values of Spatio-Temporal Entropy on 15 — 25 August, 1998 as a function of the time for times series of IMF B, , By and B- components 
in the solar wind. The thick curve represents the Entropy Index (EI) calculated over the analyzed period. We identify the MCI. EI is zero after the 
cloud due to the adjustment in the time window for calculating the entropy and that this cloud is small. 
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Figure 8: The rotation of the magnetic field vector in the plane of maximum variance and in the plane of minimum variance for the bipolar MC 
observed by ACE on 19 August at 10:00 UT to 20 August at 00:00 UT 1998. In the Figure|2]shown this MC (MCI), the EI help to identify the MC. 
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378 On 19 August 1998, at 02 : 00 UT the parameters were: proton density Np - 3.51cm , the proton temperature 

379 Tp = 7.31 * lO'* K, the ratio of alphas/protons = 4.5 * lO""*, the proton speed Vp - 333.9 km/s, magnetic field 

380 magnitude B = 2.4 nT and plasma beta /3 = 1.6. One hour later some parameters had a change to A^^ = 3.98 cm~^, 

381 Tp = 7.1 * 10"* K, alphas /protons = 6.1 * 10"^ = 331.1 km/s, B = 0.57 nT and /3 = 29.8. It is very diflicult to 

382 identify a shock, but the plasma beta increased abruptly and we believe it was related to the arrival of an event. Two 

383 hours later, at 05 : 00 UT, the plasma beta decreased to (3 - 0.58, after that there were eleven hours with gaps. The 

384 plasma beta minimum value was /? = 7.8 * 10"^ at 21 : 00 UT where there was a maximum magnetic field magnitude 

385 of B - 13.3 nT. We used the magnetic field rotation confined to one plane, the plane of maximum variance 

386 (see Figure [8]l, to find the boundaries of the cloud. 

387 Finally, we identified the bipolar MCI observed by ACE on 19 August at 10 : 00 UT to 20 August at 00 : 00 

388 UT (see Figure [8]). This MC has a flux-rope type SEN and the observed angular variation of the magnetic field is 

389 left-handed. The MVA method gives to MCI the eigenvalue ratio /I2//I3 = 18, the angle between the first and the 

390 last magnetic field vectors x - 142 °, the orientation of the axis {^Ci Oc) - (135 °, 19 °), the direction of minimum 

391 variance (0m,„, 0„„„) = (43 °, -6 °) and eigenvalues [/li, /I2, /I3] - [15.9, 12.6, 0.7]. In Table|4]we show the date of the 

392 two MCs, the time between the first and second MC is 8/1 (sheath of the second MC). The minimum speed in the first 

393 MC (MCI) was Vpi = 272.9 km/s at 15 : 00 UT and in flie beginning of the second MC (MC2) the maximum speed 

394 was Vp2 - 346.1 km/s one day later, at 05 : 00 UT. Then, Vp2 > Vp] and the two events are very close to each other. 

395 Probably there are other similar cases where a MC was confused within the sheath of the preceding MC. The entropy 
index (El) can help to identify the MCs and such kinds of cases and find the boundary of MC with MVA method. 



Table 4: Two continues MCs from 19 - 21 August 1998. The EI helps to identify the first MC. 



Year Shockl 


MCI, start 


MCI, stop 


MC2, start 


MC2, stop 


1998 19 Aug, 02:00 


19 Aug, 10:00 


20 Aug, 00:00 


20 Aug,08:00 


21 Aug, 18:00 



396 

397 3.2.3. 18-Febmary-1999 MC event 

398 A Study in more details of the STE variations is presented for some days before and after the event number 1 1 

399 Tabla[3] Time series of the magnetic components B^, By and B^ for MC have nonzero STE values (see, Figure|4^, the 

400 eleventh point). Why is this happening? To answer this question we calculate the windowed STE of three magnetic 
components B^, By and B~ in windows temporal of ten days. 




Figure 9: Values of spatio-temporal entropy at 14 - 24 February, 1999 as a function of the time for times series of IMF B^. By and B- components 
in the solar wind. The thick curve represents the Entropy Index (EI) calculated over the analyzed period. 
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Cane and RichardsonI (120031) detected one ICME not classified as MC at the date 13 Feb 19 : 00 - 14 Feb 15 : 00 
1999 with VicME = 440 km/s, Vmax - 470 km/s and B - 9 nT. Figure |9] shows the results of the calculation 
of STE from 14 Feb 15 : 00 1999 and the EI has a smafl but non-zero value. At the date 18 - 21 Feb 1999 the 
ACE spacecraft detected one ICME classified as MC with Vicme = 520 km/s, Vmax = 700 km/s, Z? = 8 nT and 
Vt - 870 km/s was the transit speed dCane and RichardsonI 2003). The MVA method gives A2/A2 - S, x - ° ^nd 
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{(pc, 6c) - (96 °, 6 °). This bipolar M C has a flux-rope type NWS and the observed angular variation of magnetic field 



is left-handed dHuttunen et al.ll2005h . Inside the MC represented by vertical lines in the Figure |9l one can see that the 
S TE = 0% only for the component. In this case, during the analyzed 10 days the EI = 0% only inside the MC. 

So far, we see MCs with zero values of STE only in one component, in two or three at a time. The interpretation 
is that it may be related to the inclination of the axis of the MC relative to the ecliptic plane and the distance from the 
spacecraft to the axis, but we have poor statistics and this issue remains open. However, the El has always detected 
the presence of a magnetic cloud. 

3.2.4. 13-October-2000MC event 

It is presented a study in more details of the STE variations with some days before and after the event number 24 
Tabla|3] Time series of the magnetic components By and B. for MC have zero STE (see. Figure]?^ the twenty-fourth 
point) but Bx has nonzero STE. 

Figure [TOl shows the results of the calculation of STE from 8-18 October 2000. On October 8 the B^, By 
components had small values of STE but diff'erent from zero. At this date it was not reported the presence of ICME, 
but we saw a small magnetic structure where the field had a smooth variation. The El detects these small structures 
that are unnoticed in most of the studies. 
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Figure 10: We show the values of spatio-temporal entropy on 08 - 18 October, 2000 as a function of the time for times series of IMF Bx, By and 
Bz components in the solar wind. The thick curve is the calculation of Entropy Index (El) over the analyzed period. 



422 In the Figure [TT] (a) solar wind parameters during one MC event are shown. From top to bottom: magnetic field 
strength, polar (Biai) and azimuthal {Bio„g) angles of the magnetic field vector in GSE coordinate system, solar wind 
speed and plasma beta. The first vertical line indicate the shock, the other two vertical lines indicate the interval of 
a MC. The plasma beta minimum value were yS = 2.7 x lO"-' on 07 October at 23 : 00 UT. We used the magnetic 
field rotation confined to one plane, the plane of maximum variance (B*B*) to find the boundaries of the cloud (see 
Table|5]l. The MVA method gives the eigenvalue ratio /I2//I3 = 47.5, the angle between the first and the last magnetic 

428 field vectors ;t' — 148.4°, the orientation of the axis (<pc,6c) - (113.8°, 7°), the direction of minimum variance 

429 (<f>„,in,Omi„) = (22°,-10°) and eigenvalues [Ai,A2,A2] = [124.6,41.8,0.9]. This MC has a flux-rope type SEN (see 
pa nel two and three in Fig urefTTIfa)') and the observed angular variation of the magnetic field is left-handed. 



ICane and RichardsonI (120030 detected one ICME classified as MC in the date 12 October 22 : 28 - 14 October 



432 17 : 00 2000 with Vicme = 4l0km/s, V,nax = 470km/s, B = 13nr and Vt = 580km/s. The MVA method gives 

433 /I2//I3 = 4, = 62 ° and (0c> 6c) = (33 °, -2 5 °). This bipolar MC h as a flux-rope type NES and the observed angular 
variation of magnetic field is right-handed ("Huttu nen et al. , 2005 1). The By, B, components have STE - 0% only 



during the passage of the MC, approximately in the first half of it and the EI detected the presence of the MC. 



Table 5: Magnetic Cloud from 07 - 08 October 2000. 

Year Shock MC, start MC, stop 

2000 07 October, 09 : 00 07 October, 22 : 00 08 October, 17 : 00 

435 
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Figure 11: Left and right part : Solar wind parameters during two MC events. Top to bottom: magnetic field strength, polar (Blat) and azimuthal 
(Blong) angles of the magnetic field vector in GSE coordinate system, solar wind speed and plasma beta. The first vertical line indicate the shock, 
the other two vertical lines indicate the interval of a MC. Left: 07 - 09 October 2000, the bipolar MC observed by ACE on 07 October at 22 : 00 
UT to 08 October at 17 : 00 UT 2000. Right: 01-11 April 2010, the bipolar MC observed by ACE on 05 April at 16 : 00 UT to 06 April at 14 : 00 
UT2010. 
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3.2.5. Interface between the interplanetary ejecta and the high-speed stream on 22 October 1999 

High-speed stream s, originating in coronal h ol es, are observed often following interplanetary coronal mass ejec- 
tion (ICME) at 1 AU JKlein and Burla^aL Il982h . IPal Lago et aP (l2QQ6h studied the 17 - 22 October (1999) solar- 
interplanetary event, which was associated to a very intense magnetic storm {Dst - -237 nT). They presented an 
analysis of pressure balance between the ICME observed on 21 - 22 October and the high-speed streams following 
this ICME. Close to the Earth, at LI, an interplanetary shock was detected by ACE magnetic field and plasma instru- 
ments on 21 October (1999), at 01 : 34 UT, as shown in Figure[T2]with the first vertical dotted line. The driver of this 
shock is an ICME, which can be distinguished from the normal solar wind b y its intense magnetic field, of the order 



of 20 nT throughout the most part of 21 October, and its low beta 0.1) ( Dal Lago et al. . 20061). The start of the 



ejecta was at 03 : 58 UT of 21 October. Toward the end of this ejecta, an increase of the magnetic field intensity was 



observed, starting at 02 : 30 UT of 22 October, reaching a peak value of 37 nT dDal Lago et al. , l2006l) . At 06 : 15 



UT of 22 October (th e second vertical dotted line in Figure [T2]). the magnetic field dropped abruptly around 10 nT. 
Dal Lago et al. (l2006h defined this p oint as the end of the ICME. They were not sure whether this ICME is an MC 



or not acc ording to the crit eria of Burla ga et al.l (119811) . because the direction of the magnetic field does not rotate 
smoothly. ICane and Richardson (.2003.) also detected one ICME not classified as MC at the date 21 October 08 : 00 
UT - 22 October 07 : 00 UT 1999 with Vjcme = 500 km/s, V,„„^ = 580 km/s, B = 20 nT and Vj = 480 Km/s. 

As presented in the previous paragraph it is clear that during the period 20 - 26 October 1999 was not identified 
any magnetic cloud. In Figure [121 we show the values of spatio-temporal entropy as a function of the time for times 
series of IMF B^, By and B. components in the solar wind from 20 - 26 October. All time the values of STE are 
different of zero, so the El is also different from zero. This result helped to validate the El to detect MC, because El 
methodology can differentiate when an ICME is not classified as a MC, since the El is zero only inside a MC. 



100 


« — 


-ICME > 






ACE 


80 
60 










Bz4 ; 


40 










20 























GSM 






Oct 21 


Oct 22 


Oct 23 


Oct 24 


Oct 25 



1999 



Figure 12: Values of spatio-temporal entropy on 20 - 26 October, 1999 as a function of the time for times series of IMF B, , By and B- components 
in the solar wind. The thick curve represents the Entropy Index (EI) calculated over the period analyzed. 



457 3.2.6. Identification by windowed STE method of 05-06 April 2010 MC event 

458 On 3 April 2010 the Sun launched a cloud of material, known as a coronal mass ejection (CME), in a direction 

459 that reaches the Earth. This CME was much faster than most of the CMEs seen recently, with a speed of at least 

460 800 Km/s. The bulk of the CME passed south of Earth, but a piece of it hit the Earth's magnetosphere on April 5, 

461 causing a geomagnetic storm {Dst^y^ - -73 nT on 06 April at 15 : 00 UT). As a result bright auroras were seen 

462 around the world at higher latitudes. 

463 This event was very recent and we wanted to use the technique to verify any possible MC occurrence. The ACE 

464 Magnetic Field Experiment data in level 2 (verified) were not available in April 2010 when the data were processed. 

465 It was only possible to obtain such data for 16 second average IMF in RTN and GSE coordinates via anonymous ftp, 

466 from Caltechlfl For this reason, in Figure[T3]the STE values are calculated for data in GSE coordinates. 



^ftp://mussel. srl.caltech.edu /pub/ace/browse/MAG16sec 
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Figure 13: Values of spatio-temporal entropy on 01 - 11 April, 2010 as a function of the time for times series of IMF B^. and B, components 
in the solar wind. The thick curve represents the Entropy Index (EI) calculated over the period analyzed. The MC is identified. 



FigurefTSlis similar to Figure|6]but for the period 01 - 10 April 2010 with GSE coordinates. We find minimum 
values of S TE - l%to B, on 05 April at 21 : 33 : 20 and 22 : 26 : 40 UT respectively. The and By components 
had STE values less than 10% on 05 April at 18 : 00 : 00 and 18 : 53 : 20 UT respectively. Then, the EI has a 
minimum values less than 1% on day 05 between 18 : 00 : 00 to 23 : 20 : 00 UT. Thus the EI de tects a structure with 
characteristics of a magnetic cloud. The same method (MVA) explained in lHuttunen et al.l(l2005h was used to find the 
boundaries of the cloud. 



Table 6: Magnetic Cloud from 05 - 06 April 2010. 

Year Shock MC, start MC, stop 

2010 05 April, 07 : 00 05 April, 16 : 00 06 April, 14 : OO" 

472 

473 Figure [IHb) is similar- to (a) but from 01-11 April 2010. On 05 April 2010, at 07 : 00 UT, the proton density 

474 Np - 2.8 cm""*, proton temperature Tp - 2.4 * 10^ K, ratio of alphas I protons = 7 x 10~^, proton speed Vp - 

475 564.5 km/s, magnetic field magnitude B - 5.2 nT and plasma beta fi — 0.9. One hour later some parameters have an 

476 change = 7.98 cm-^ Tp = 5.5 x 10^ K, alphas I protons = 1.7 * 10"^ Vp = 724.8 km/s, B = 10.9 nT and/? = 1.3. 

477 It is easy to identify a shock because the velocity, density and magnetic field magnitude increases abruptly and we 

478 believe it is related to the arrival of an event to ACE. Eight hours latter than the time of the shock, at 16 : 00 UT, the 

479 plasma beta decreases to - 3.2 * 10"^, and it is the beginning of the magnetic cloud. The plasma beta minimum 

480 value were [i - 8.8 x 10"^ on 06 April at 12 : 00 UT. We used the magnetic field rotation confined to one plane, 

481 the plane of maximum variance {B*B*) to find the boundaries of the cloud (see Table|6]). The MVA method gives the 

482 eigenvalue ratio /I2//I3 = 16, the angle between the first and the last magnetic field vectors^ - 48 °, the orientation 

483 of the axis {(pc,0c) - (33 °, -32.6 °), the direction of minimum variance {<pmi„,Omin) = (129°, 9°) and eigenvalues 

484 [Ai,A2, A3] - [16.8, 3.6, 0.2]. This MC has a flux-rope type NWS (see panel two and three in Figure [TT](b)) and the 

485 observed angular variation of the magnetic field is left-handed. 

486 3.3. Correlation between the Windowed STE and the inclination of IMF 

487 Following the same idea presented in subsection l2.3l we calculate the STE versus the trend angle. In Figure[3] we 

488 presented time series plot of Lorenz data file include in the VRA for four differents trend angles. We mentioned earlier 

489 that the time series of magnetic clouds are not stationary, there is a trend. The linear least squares fitting technique is 

490 the simplest and most commonly applied form of linear regression and provides a solution to the problem of finding 

491 the best fitting straight line through a set of points and to obtain the equation of the line. The angle 6 that a line makes 

492 with the positive x axis is closely related to the slope m via the inverse tangent function 9 - tan"'(m). We applied the 

493 previous methodology to all IMF time series shown in the Figures |6l|7]|9l[l0l[l2l In other word, we separate all pairs 

494 of points formed by the ti^end angle and STE for the data. 
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Figure 14: At top: Plot of STE versus 6 to 6,-, 6,. and B. components and a parabola fitted to the data. At bottom: The corresponding histogram in 
which a Gaussian distribution was fitted. 
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495 Figure [T4lpresents for the Bx, By and Bz components, at the top of each panel, the plot of STE versus 6 for each 

496 magnetic component and a parabola fitted to the data, and at the bottom of each panel, the corresponding histogram, in 

497 which a Gaussian distribution was fitted. At the first panel in the figure, for Bx, from the STE versus 0, the correlation 

498 coefficient (a Pearson's correlation coefficient between two time series, STE vs 6) found was C - -70.8%. If the 

499 variables tend to go up and down in opposition, with low values of one variable associated with high values of the 

500 other, the correlation coefficient will be negative (anti-correlation), as in our cases. The anti-correlation between the 

501 two parameters (STE and 6) gives the idea that the trend in the Bx series is responsible for the low values of STE. But 

502 the above idea is not totally true; because we have S TE + 0% series with a larger trend, more separated points of the 

503 vertical line. The histograms help to see the distribution of the tend. 

504 One response to these kinds of problems can be fluctuations in the signal. Sometimes, the trend in the records for 

505 a sheath region is bigger than the trend for a cloud, but the sheath has rapid fluctuations that changes the sign of the 

506 signal and the STE is not zero. Then the S TE - 0% at magnetic clouds is related to the simultaneous occurrence 

507 of two factors: (1) the trend existing in the signal, (2) the analyzed variable (in this case, IMF B^) has a more ordered 

508 dynamical behavior {i.e. few fluctuations) and higher degree of correlation between its time neighbors. 

509 In Figure[T4]the middle panel, for By, and the bottom panel, for Bz, present respectively the correlation coefficient 

510 (C) results of -67.6% and -65.7%. From the three panel, it is possible to see that B^ component has largest trend 

511 compared to the others two, but the reason is not know. 

512 4. Conclusions 

513 We dealt with time series of solar wind for a group of magnetic clouds in order to analyze the dynamical behavior 

514 of the IMF Bx, By and B- components. The tool chosen to this study was the spatio-temporal entropy (called STE) 

515 method. From methodological considerations based on physical reasons, the STE was calculated from solar wind time 

516 windows with 2500 points corresponding to 11.11 hours of data. The results were significant and very promissing. we 

517 found higher entropy values in the sheaths. MCs present zero STE values at least for one magnetic component. As not 

518 all the magnetic components have zero entropy at the same time, an useful technique was to create a standardization 

519 index, called Entropy Index (El), to allow joining the STE results of the three variable (B i, By, B.). We studied in 

520 more detail the STE variations taking into account some days before and after to six selected events, in order to allow 

521 understanding and validating the analysis procedure. By this approach, the El identified three new MCs in addition to 

522 the known MCs presented in previous tables from other researchers, and after the MVA method was used to calculate 

523 the boundaries easily. 

524 The analyses developed in this work shows that: 

525 • MCs have STE values lower than the ones for sheath region. 

525 • The differences among the STE values for the three magnetic components in a MC give an idea about the 

527 anisotropy in the structure of some MCs. But this idea is still an open question to be better investigated. 

528 • In MC the magnetic field strength is higher than the average, the magnetic field direction rotates smoothly 

529 through a large angle, then periods of MC have more trend than sheath region and "quiet" solar wind periods. 

530 The trend is the principal cause of the lower values of STE in the MC. Also, MCs are most structuredness than 

531 sheath and "quiet" solar wind. 

532 • Using the EI and MVA method it is possible to identify all the magnetic clouds that arrive at spacecrafts. 

533 • The entropy index (EI) proposed gets success as an auxiliary tool to identify magnetic cloud candidates and 

534 calculation of EI could be used in an automatic procedure of preliminary investigation or convenient warning. 

535 Based on its skill to provide an objective analysis, the STE method can be seen as a valuable tool in the study of 

536 interplanetary phenomena such as the time localization of MCs. 
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